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ABSTRACT 
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^ , Recent investigations of the magnetic field vector properties in the solar internetwork have provided diverging results. 

. While some works found that the internetwork is mostly pervaded by horizontal magnetic fields, other works argued 

■ in favor of an isotropic distribution of the magnetic field vector. Motivated by these seemingly contradictory results 
and by the fact that most of these works have employed spectropolarimetric data at disk center only, we have revisited 
this problem employing high-quality data (noise level cr ~ 3 x 10~* in units of the quiet-Sun intensity) at different 
latitudes recorded with the Hinode/SP instrument. Instead of applying traditional inversion codes of the radiative 

■ transfer equation to retrieve the magnetic field vector at each spatial point on the solar surface and studying the 

■ resulting distribution of the magnetic field vector, we surmised a theoretical distribution function of the magnetic field 
vector and used it to obtain the theoretical histograms of the Stokes profiles. These histograms were then compared 

(~| ■ to the observed ones. Any mismatch between them was ascribed to the theoretical distribution of the magnetic field 

^ ^ vector, which was subsequently modified to produce a better fit to the observed histograms. With this method we find 

I ■ that Stokes profiles with signals above 2 x 10~^ (in units of the continuum intensity) cannot be explained by an isotropic 

O I distribution of the magnetic field vector. We also find that the differences between the histograms of the Stokes profiles 

observed at different latitudes cannot be explained in terms of line-of-sight effects. However, they can be explained by a 



CZ2 ■ distribution of the magnetic field vector that inherently varies with latitude. We note that these results are based on a 



series of assumptions that, although briefly discussed in this paper, need to be considered in more detail in the future. 
Key words. Magnetic fields - Sun: photosphere - Sun: surface magnetism - Sun: magnetic topology 
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QQ . 1. Introduction of the magnetic field vector with respect to the observer's 

■ hne-of-sight), others favor a quasi- isotropic distribution of 

Q ■ In recent years several attempts have been made to magnetic fields (Martinez Gonzalez et al. 2008; Asensio 

• ; mvestigate the magnetic field vector distribution m the j^^mos 2009; Stenflo 2010). With a few exceptions (Harvey 

^ . solar internetwork. Initially, these works studied the al. 2007, Lites et al. 2008 and Martinez Gonzalez et al. 

^ ■ magnetic field strength is m these regions. Some favored 2OO8), all previous studies were carried out employing data 

magnetic fields of about a few hundred Gauss or less recorded at disk center only. Therefore, to better constrain 

(Asensio Ramos et al. 2007; Lopez Ariste et al. 2007; t^^g angular distribution of the magnetic field vector in 

Orozco Suarez et al. 2007a, Orozco Suarez & Bellot internetwork, we considered spectropolarimetric data 

Rubio 2012) while others found magnetic fields in the recorded at different positions on the solar disk (Section EJ. 
kilo-Gauss range (Dominguez Cerdeiia et al. 2003, 2006; 
I Sanchez Almeida 2005). These studies were carried out 

mostly with low spatial resolution data (1"). Whenever In addition, Asensio Ramos (2009), Stenflo (2010), 

the spatial resolution increased to better than 1 arcsec, and Borrero & Kobel (2011; hereafter referred to as 

this decreased the signal-to-noise ratio. With the Hinode paper I) warned that the highly inclined magnetic fields 

satellite (Kosugi et al. 2007) it is now possible to obtain obtained by some studies could be caused by the noise in 

spectropolarimetric data (full Stokes vector) with high the linear polarization profiles. This yields a distribution 

spatial resolution (0.3") and low noise (cr ~ 10"'^ in units of B± (component of the magnetic field vector that is 

of the continuum intensity). Thanks to these new data, perpendicular to the observer's line-of-sight) with a peak at 

it is now also possible to investigate not only the module around 50-90 Gauss. To avoid this problem, these authors 

but the three components of the magnetic field vector, proposed to include only those profiles in the analysis that 

This has led to a new controversy about the angular have a signal-to-noise ratio > 4.5 in the linear polarization 

distribution of the magnetic field vector in the quiet Sun. (Stokes Q and U). Although this selection criterion allows 

While some authors (Orozco Suarez et al. 2007a, 2007b; one to retrieve reliable distributions for the magnetic field 

Lites et al. 2007, 2008) found that the magnetic field is vector, it has the disadvantage of excluding most of the 

mostly horizontal (7 90°; with 7 being the inclination Stokes profiles within the field-of-view from the analysis 
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(see Borrero & Kobel 2012; hereafter referred to as paper 
II; cf. Bellot Rubio & Orozco Suarez 2012). In this paper 
we adopt an ahernative approach based on inverting the 
histograms of the observed Stokes vector (Section 2]) over 
the entire field-of-view instead of inverting the Stokes vec- 
tor at each pixel over the observed region. Under a number 
of simplifying assumptions, whose limitations are described 
in Section [21 we were able to reach some important, albeit 
preliminary, conclusions about the angular distribution of 
the magnetic field vector in the solar internetwork and its 
variation across the solar disk (Section [5]) . 



2. Observations and datasets 

The data employed in this work were recorded with the 
spectropolarimeter (SP; Ichimoto et al. 2008) attached to 
the Solar Optical Telescope (SOT; Tsuneta et al. 2008, 
Suematsu et al. 2008, Shimuzu et al. 2008) onboard the 
Japanese spacecraft Hinode (Kosugi et al. 2007). The 
spectropolarimetric data comprise the full Stokes vector 
(/, Q,U,V) around the pair of magnetically sensitive 
spectral lines Fel 6301.5 A (c/off = 1-67) and Fel 6302.5 
A (5eff = 2.5). QoS refers to the effective Lande factor 
calculated under LS coupling. The spectral resolution 
of these observations is about 21.5 mA per pixel, with 
112 pixels in the spectral direction. The spatial resolu- 
tion of Hinode/SP observations is 0.32". For this paper 
we selected three maps at three different heliocentric 
positions. In all three maps the spectrograph's slit was 
kept at the same location on the solar surface for the 
whole duration of the scan. This means that, while the 
vertical direction (F-axis or direction along the slit) 
contains information about different spatial structures 
on the solar surface, the horizontal direction (X-axis or 
direction perpendicular to the spectrograph's slit) samples 
the same position at different times. Each spectrum was 
recorded with a 9.6 seconds exposure, yielding a noise of 
about a = 7.5 x 10~^ in units of the quiet-Sun continuum 
intensity. Each map records data for a period of time (> 1 
hr) that includes several turnovers of the granulation, 
thus breaking down the temporal coherence and providing 
spatial information (in a statistical sense) along the X-axis. 

In paper I we have demonstrated that photon noise 
plays an important role in determining the magnetic field 
vector from spectropolarimetric observations. To further 
decrease the level of noise in our observations we aver- 
aged every seven slit positions (temporal average of about 
67.1 seconds), which yields a new noise level of about 
(T = 3 X lO^"' (in units of the quiet-Sun continuum inten- 
sity). However, averaging means that the original map is 
shortened by a factor of seven in the direction that is per- 
pendicular to the slit (X-axis). This decreases the number 
of points available for statistics. Fortunately, Hinode/SP 
data have a sufficient number of pixels to ensure good statis- 
tics even after averaging (see Section In the following we 
briefiy describe each map individually. 

2.1. Map A 

This map was recorded on February 27, 2007 between 00:20 
UT and 02:20 UT. It originally consists of 727 slits posi- 



tions, of which 103 remain after temporal averaging. The 
center of slit was located at approximately the following co- 
ordinates on the solar surface: X = —31.7" and Y — 7.7". 
This corresponds to a heliocentric position of /i = cos O ss 1 
(O is the heliocentric angle) and to a latitude of A « 0°. 
The noise level is <t = 2.8 x 10^"*. This map (original and 
temporally averaged) corresponds to Maps B and C in pa- 
per I, and it was also employed (with and without temporal 
averaging) by Lites et al. (2008) and Orozco Suarez et al. 
(2007a). 

2.2. Map B 

This map was recorded on February 6, 2007 between 
11:33 UT and 15:51 UT. It originally consists of 1545 slits 
positions, of which 222 remain after temporal averaging. 
The center of slit was located at approximately the 
following coordinates on the solar surface: X = 493.6" and 
Y = 491.3". This corresponds to a heliocentric position of 
fi — cosO « 0.7 and to a latitude of A « 30°. The noise 
in this map is very similar to that in Map B: cr = 3.1 x 10~^. 



2.3. Map C 

This map was recorded on January 17, 2007 between 
07:05 UT and 09:58 UT. It originally consists of 1048 slits 
positions, of which 149 remain after temporally averaging. 
The center of slit was located at approximately the 
following coordinates on the solar surface: X = —3.0" and 
Y — 697.1". This corresponds to a heliocentric position of 
/i — cos 6 sa 0.7 and to a latitude of A « 40°. Here, the 
noise level is slightly higher than in Map A: cr = 3.2 x 10"'*. 
We note that some consecutive slit positions in this map 
show very high noise in the Stokes profiles. Although 
we could not relate this effect to the South Atlantic 
Anomaly (increased flux of cosmic rays at certain orbits 
of the satellite) we have removed these slit positions from 
our analysis, which reduced the effective number of slit 
positions to 120. 

Figures [U [21 and [3] display the magnetic flux density 
(f> of maps A, B, and C as obtained through the inversion 
of the full Stokes vector with the VFISV (Very Fast 
Inversion of the Stokes Vector) inversion code (Borrero et 
al. 2010). For better visualization the maps in these figures 
are obtained from the inversion of the original data (i.e. 
not temporally averaged). This avoids pixelization in the 
X-axis of these plots. However, for the remainder of the 
paper, our discussions and figures are based only on the 
temporally averaged (67.1 seconds) data. 

Although these previous figures only show the to- 
tal magnetic flux density, it is worth mentioning that 
the VFISV code also retrieves the three components 
of magnetic field vector B: B is the module of B, 7 
is the inclination of B with respect to the observer's 
line-of-sight, and is the azimuth of B in the plane that is 
perpendicular to the observer's line-of-sight. In addition, 
VFISV retrieves the magnetic filling factor a as well as 
the line-of-sight component of the velocity vector Vios and 
a set of thermodynamic parameters T. We note that the 
magnetic flux density is deflned as <& = aB. For a more 
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Fig. 1. Inferred magnetic flux density $ ~ aB from the 
inversion of Map A (Sect. I2.ip . White areas correspond to 
regions where ah three polarization profiles (Stokes Q, U, 
and V) are below the 4.5tT-level. 



Fig. 3. Same as Figured] but for map C (Sect. 



Magnetic flux density [Mx cm 



150 300 450 

Map B: A = 30° ; yu. 



600 

= 0.7 

r 



750 



550 - 



^ 500 - 



fi 450 - 



■a 
I 




100 150 
Time [min] 



Fig. 2. Same as Figure [T] but for map B (Sect. 12.2 



3. Stokes profiles at different positions on the solar 
disk. 

The inversions carried out in the previous section could 
be employed to obtain histograms of the magnetic flux 
density $, module of the magnetic field vector B, and 
the inclination of the magnetic field vector with respect 
to the observer's line-of-sight (7) at different positions on 
the solar disk. However, in this paper we aimed to infer 
properties about the distribution of the magnetic field 
vector by directly studying the histograms of the Stokes 
profiles. Figure 0^ presents distribution histograms of 
the maximum signals of the Stokes V{X) (dashed lines) 
and Stokes (5(A) and U{X) (solid lines) normalized to the 
— average quiet- Sun intensity over the entire map: Iqs- The 
colors indicate each of the different maps studied: red for 
map A (Sect. 12. ip . green for map B (Sect. 12. 2p . and blue 
for map C (Sect. 12. 3p . Figure SJd displays the cumulative 
histogram of the pixels in each map that have a S/R 
(signal-to-noise ratio) equal to or higher than a given 
value. The colors and the line-styles are as in Figure |3^. 
For instance: 31.6 % of the pixels in map A posses signals 
in (5 or [/ (solid-red line) that are above 4.5 times the noise 
level. To limit our analysis to the internetwork regions, we 
excluded from these figures the pixels in maps A, B, and 
C with a magnetic flux density $ > 500 Mx cm~^. 



250 



detailed overview on Milne-Eddington inversion codes, 
which include not only the magnetic field vector but also 
the thermodynamic and kinematic parameters relevant to 
the line formation, we refer the reader to del Toro Iniesta 
(2003), Borrero et al. (2010) and references therein. 



We now focus on some of the features of the histograms 
in Figure m A very striking one is the peak at around the 
3(7- value for the linear polarization (solid lines) in Fig. 0^, 
and the fact that the amount of pixels with maximum linear 
polarization signals below and above this 3(T-level quickly 
drops. One might wonder how a peak can appear at around 
3(7 if the probability that photon noise (taken as a normal 
random distribution) will produce such a high value is only 
about 0.27 % ? The answer to this question is that photon 
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Fig. 4. Left panel: histogram of the number of pixels as a function of the maximum value in their polarization signals 
(normalized to the average continuum intensity). The vertical solid-black line represents the Scr- level. Right panel: his- 
togram of pixels with polarization signals above a certain signal-to- noise ratio (S/R). This last histogram is cumulative, 
while the first one is not. In both panels the color lines indicate the same: red (map A at A = 0°), green (map B at 
A = 30°), blue (map C at A = 40°). The solid lines refer to the linear polarization profiles (Stokes Q and U) while 
the dashed lines correspond to the circular polarization profiles (Stokes V). For comparison we also display, in solid and 
dashed-black curves, the expected behavior due to noise. 



noise has a probability of 0.27 % to produce a signal at the 
3(T-level at one particular wavelength, but Figure 2] shows 
the maximum of the signal over all wavelengths. Indeed, 
it is possible to employ the binomial distribution to find 
the lower bound of the probability K{d, N) that after N 
wavelengths, one of them will have a signal stronger than 
or equal to 5-times the noise level a: 

K{6,N)>l-[l-p{6)]'' , (1) 

where p{S) is the probability that at one single wavelength 
position, a normally distributed random variable (with a 
standard deviation a) will yield a signal (5-times above the 
standard deviation is given by 

p{S) = 1 - erf(5/V2) , (2) 

where erf denotes the so-called error function. Since the 
spectral line is sampled in TV = 112 spectral positions for 
each Stokes parameter, the probability of finding a wave- 
length where the noise yields a signal at the 3ct is much 
higher than the 0.27 % mentioned above. In particular, 
there is a K{N = 112,(5 = 3) > 26.12% probability that 
the noise in the circular polarization (Stokes V) will yield a 
signal at the 3cr-level. Because the linear polarization con- 
sists of two Stokes parameters {Q and U), this probability 
is even higher: K{N = 224, ^ = 3) > 45.42%. Given these 
high probabilities, it is not surprising that the histograms 
of the Stokes profiles peak close to the 3o'-level. This is 
certainly the case of the linear polarization (solid lines in 
Fig. H^). It is noteworthy that, at the 3cr-location, Stokes 
V (circular polarization) only presents a local maximum 
(dashed lines in Figure H^) that actually disappears for 
maps B (dashed green) and C (dashed red). Indeed, the 
peak for the histogram of Stokes V appears to be located 
in the range of 2 — 3 x 10"'^ (in units of the quiet Sun 
continuum intensity), which is about 8 times above the 



noise level a. 

These differences between linear and circular polar- 
ization in Figure 0^ can be explained if we consider 
a distribution of (component of the magnetic field 
vector that is parallel to the observer's line-of-sight) that 
features a peak at a value such that the corresponding 
Stokes V profile would have a maximum at around Scr. 
At the same time, the distribution of B± (component 
of the magnetic field vector that is perpendicular to the 
observer's line-of-sight) should feature a low probability 
of finding values of B± that produce Stokes Q and U 
profiles above 3a such that there is a peak at this level. 
We additionally investigated this by synthesizing Stokes 
profiles with different values of -B|| and B± and comparing 
them with the maximum signal of the resulting Stokes 
profiles. To obtain a good estimation we carried out this 
experiment with two different semi-empirical models that 
represent granules and intergranules (Borrero & Bellot 
Rubio 2002). Results are presented in Figure [SJ This figure 
shows that the distribution of that is responsible for the 
Stokes V signal in Figure (dashed-lines) must posses a 
peak at around -B|| ss 5 — 7 G, which is the value needed to 
produce a majority of Stokes V signals at the 8cr-level. In 
addition, Figure [S] shows that the distribution of B± must 
have very low probabilities for B± > 40 — 70 G, otherwise 
the peak in the histograms for Q and U in Figure (solid 
lines) would be significantly shifted above the 3cr-level. 

An important feature of Figures and |3)d is that the 
distribution of the linear polarization signals (solid lines) 
is different at different positions on the solar surface. In 
addition, the distribution of the circular polarization sig- 
nals (dashed lines) also changes slightly, but comparatively 
less than the linear polarization. 
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Fig. 5. Amount of circular polarization (normalized to the 
quiet-Sun continuum intensity) generated as a function of 
the component of the magnetic field vector that is parallel 
to the observer's line-of-sight: (dashed lines). Amount of 
linear polarization (normalized to the quiet-Sun continuum 
intensity) signals generated as a function of the compo- 
nent of the magnetic field vector that is perpendicular to 
the observer's line-of-sight: i?_L (solid lines). These curves 
were obtained by performing a synthesis of the Stokes pro- 
files employing different values of and and using 
two different atmospheric semi-empirical models: red for 
intergranules and blue for granules. The horizontal black 
lines indicate the 2>a and Scr levels assuming a noise of 
cr = 3 X 10~^ as in our observed maps (Section . 



4. Results from theoretical distributions 

In the previous section we have determined some general 
properties about the distribution of the magnetic field 
vector in the internetwork by considering some details 
from the observed histograms of the Stokes profiles. We 
now continue along these lines by investigating the sources 
of the differences in the distribution of the polarization 
signals at different positions on the solar disk. Following 
the notation introduced in paper II, we refer to X as the 
set of physical parameters that affect the formation of the 
Stokes profiles: 



X= [T,yios,B,a] 



(3) 



where T refers to the thermodynamic parameters, T^os to 
the line-of-sight-velocity, B to the magnetic field vector, 
and finally a refers to the magnetic filling factor (see also 
Section [2]). Considering a Milne-Eddington atmosphere 
(see del Toro Iniesta 2003; Borrero et al. 2010) is equivalent 
to assuming that with the exception of source function, 
none of the thermodynamic, kinematic, and magnetic 
parameters change with optical depth in the photo- 
sphere. The source function, however, is considered to 
vary linearly with optical depth: S(t^ = 5*0 + t^Sx, where 
5*0 corresponds to the source function at the observer's 
location, and S\ is the derivative of the source function 
with optical depth. In addition to and 5*1, the other 
thermodynamic parameters in T are the Doppler width 



of the spectral line A Ad, the damping parameter a, and 
quotient of the absorption coefficient in the continuum and 
in the line center 770- 

We now refer to 7^(X)dX as the probability of finding 
a pixel within the observed field-of-view where each of 
the physical parameters in X has values between Xi and 
Xi -\- dXi. Furthermore, we assume that the magnetic 
parameters are statistically independent of the thermody- 
namic and kinematic parameters, thereby allowing us to 
write 



7'(X)dX = 751(3, 0)7^2 (T,Mos)dBdTdadl1os • 



(4) 



We now turn our attention to the probability distribu- 
tion function of the magnetic parameters Vx , which can be 
rewritten as 



751(6, a)dBda =ri{B„ Bj, Bk, a)dB,dBjdBkda , 



(5) 



which indicates the probability of finding a pixel whose 
magnetic field vector has the coordinates between the fol- 
lowing values: Bi and Bi + dBi, Bj and Bj + dBj, Bk and 
Bk + dBk, and finally where the filling factor of the mag- 
netic field has a value between a and a + da. The reason 
for our choice of nomenclature in the three components of 
the magnetic field vector Bi, Bj and Bk is that the proba- 
bility distribution function will generally not be expressed 
in spherical coordinates in the observer's reference frame 
(the frame needed to solve the radiative transfer equation) . 
In general, we therefore must perform a transformation of 
variables that will express Equation [S] into the observer's 
reference frame and into spherical coordinates: 



ri{B„Bj,Bk,a)dB,dBjdBkda^\J\Vi{B,-f,(j),a)- 

■ dBd'yd(l)da , 



(6) 



where |J| is the determinant of the Jacobian matrix for 
the transformation between the two reference frames. This 
transformation might also introduce a dependence on the 
heliocentric angle Q, thereby allowing us to evaluate the 
probability distribution function at different positions on 
the solar disk. The total probability must be equal to one: 



da 



dB 



27r 



d7 / |J|7'i(i?,7,</',a)d</)=l . (7) 



We then use this distribution function to evaluate the 
percentage of pixels (from the total) that posses a given 
magnetic field vector. In our simulations we employed 
a total of 2 X 10^ pixels. To solve the radiative transfer 
equation, we need the probability distribution function 
of the thermodynamic and kinematic parameters V2 (see 
Eq. |4]) , in addition to the probability distribution function 
of the magnetic parameters Vi- Hereafter we take, as 
7^2(T, Vios)dTdVios, the distribution obtained from the 
results of map A (called map C in paper I) . We assume that 
this distribution does not depend with the position on the 
solar disk. Finally, in all our tests in this section we assume 
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that the magnetic filhng factor is unity: a — 1. We take 
this approach not to add a new degree of freedom that will 
make our subsequent analysis more cumbersome. Section [5] 
will address all the assumptions and simplifications of this 
section in more detail. 

With this we now have all the necessary ingredients 
needed by the VFISV (Borrero et al. 2010) code to solve the 
radiative transfer equation in order to obtain theoretical 
Stokes profiles Q, U, and V. To these profiles we then add 
noise, assuming a normally distributed random variable 
(Leva 1992) with a standard deviation cr = 3 x 10^* (in 
units of the quiet-Sun continuum intensity). This noise 
level is similar to that of the observed maps in Section (3] 
Once noise has been added, we select the peak values of the 
Stokes profiles and construct histograms like those derived 
from the observations in Figs. |3^-|3]3. The theoretical and 
observed histograms are then compared with different 
theoretical distributions of the magnetic field in a attempt 
to match the observations. 

The approach described in this section is indeed an 
inversion. However, it is not the same kind of inversion as 
those in Section[2]or in paper I. First of all, the observables 
here are the histograms of the peak values in Stokes Q, U, 
and V (Figs. |4^-|4)d), whereas in paper I the observables 
were the full Stokes vector (/, Q, U, and V) including 
their wavelength dependence (not only the peak values) at 
each individual pixel. The model parameters also differ: 
while before the model parameters corresponded to the 
three components of the magnetic field vector for each 
individual pixel, here the model parameters correspond 
to a parametrized theoretical distribution function for the 
magnetic field vector that includes all pixels. As we show 
below, this parametrized theoretical distribution function 
possesses a very limited number of free parameters (1-3), 
which can be tuned to simultaneously fit all pixels in the 
field-of-view and at different heliocentric angles. This is 
a great advantage compared to traditional inversion of 
Stokes profiles (e.g. Borrero & Kobel 2011; paper I), where 
there were about ten free parameters (those describing 
a Milne-Eddington atmosphere) for each pixel in the 
map. Another important difference is that, while regular 
inversion codes for the radiative transfer equation (Ruiz 
Cobo & del Toro Iniesta 1992; Borrero et al. 2010) are 
automatized, the procedure followed here is completely 
manual (i.e: trial and error). 



4.1. Isotropic distribution functions 

In this section we employ a theoretical distribution where 
the magnetic field vector is isotropic. We first define as 
isotropic a probability distribution function where the mag- 
netic field vector has no preferred orientation. In the local 
reference frame on the solar surface this can be expressed 
as having a probability of finding a magnetic field vector 



that depends only on its module B = \/B^ 



Bh 



Pi(B)dB = Af{B)dBABydB^ 



(8) 



where A is just a normalization constant. In this reference 
frame, {ex, By, e^}, the ez-axis is perpendicular to the solar 
surface. To take into account that the observer's line-of- 



sight forms an angle Q (heliocentric angle) with respect 
to the Gz-axis, we perform a variable change into a new 
coordinate system {ex,ey,ez} that is rotated by an angle 
Q around the ey-axis. Indeed, because the distribution is 
isotropic, it does not matter around which axis we consider 
the rotation. With this transformation the new e^-axis is 
aligned with the observer's line-of-sight. In this case, the 
relationship between the old coordinates of the magnetic 
field vector with the new ones is given by 



B^ = sin e + cos e , 
B2 = B^ cos Q — B^smQ , 

Ul = 1 . 



(9) 



Since this is a simple rotation, the determinant of the 
Jacobian matrix for the transformation is unity. In addi- 
tion, the module of the magnetic field vector is the same in 

the old and new reference frames: B = ^/B'^ + B'^ + B^ = 



B'^ + B'^ -I- B^^. We can therefore rewrite Equation [5] in 
the observer's reference frame as 

Pi(B)dB = A/(B)dB^dBydB^ , (10) 

which is functionally identical to Eq. |S] since it was defined 
to be isotropic and therefore independent of the viewing 
angle 0. As previously mentioned, to solve the radiative 
transfer equation we need to express Equation [TU] in spher- 
ical coordinates. We therefore perform now an additional 
variable change where 



B^ = B sin 7 cos < 
By ^ B sin 7 sin (, 
B^ — B cos 7 , 
|J| = 52 sin7 . 



(11) 



and thus. 



Pi(B)dB = Af[B)B'^ sin7dBd7d0 



(12) 



Since the distribution is isotropic, it must also be inde- 
pendent of Q once it is expressed in the observer's refer- 
ence frame both in Cartesian coordinates or spherical ones 
fEas.[T0land[T2l respectively). In addition, isotropism man- 
ifests itself as a sin7-dependence, which comes from the 
determinant of the Jacobian matrix or, in order words, the 
volume-element expressed in spherical coordinates: dB = 
B^ sin 7di3d7d0. For the distribution of the module of the 
magnetic field we consider an exponential function in the 
form f{B) « exp(— _B). After normalization, the resulting 
expression for the theoretical distribution function is 



'Pi{B,-i,4>)dBA-id(l)= ^;^exp 



27 B2 

'0 



-3B 



sin 7d-Bd7d(/) , 
(13) 

where Bq represents the mean value of the magnetic field 
vector module: 



=< B > = 



Jo 



2ir 



B ■Vi{B,'y,(t))dBd-fd(j) . (14) 
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Fig. 6. Panel a: histogram of the number of pixels with a given maximum value in their polarization signals (normalized 
to the average continuum intensity in the quiet Sun). The vertical solid-black line corresponds to the 3cr-level. Panel b: 
histogram of pixels with polarization signals above a certain signal-to-noise ratio (S/R). This last histogram is cumulative, 
but the first one is not. The solid and dashed lines refer to the linear and circular polarization profiles, respectively. These 
two panels can be readily compared to the observed histograms in Figure [H They were obtained employing an isotropic 
distribution function for the magnetic field vector (Eq. [T3|) with three different mean values for the magnetic field strength 
Bq: 20 (red), 30 (blue), and 40 (green) Gauss. The module of the magnetic field vector B and the inclination of this 
vector with respect to the observer's line-of-sight 7 arising from these isotropic distributions are presented in panels c 
and d. 



Figures Ht-IBJl display the probability distribution 
functions for module B and inclination of the magnetic 
field vector with respect to the observer's line-of-sight 7, 
corresponding to Eq. I13[ and employing three different 
values of Bq — 20, 30, and 40 G. The histograms for 
the resulting Stokes profiles after solving the radiative 
transfer equation with these distributions are given in 
Figures in^-O), which can be readily compared to Figs. |3^- 
[4)d. This comparison shows that for Bq ~ 20 G, the linear 
polarization produced by an isotropic distribution (solid 
lines) is too low compared to the observed one. For Bq = 40 
G, the amount of linear polarization is comparable to the 
observed one, but, in this case the isotropic distribution 
produces too much circular polarization (dashed lines). In 
addition, an isotropic distribution produces far too few Q 
and U profiles, about 50 % of the total, with S/R > 3 (sohd 
lines in Fig. |B]d), while the observed distribution shows 
that about 70-80 % of the profiles are above tliis level 



(solid lines in Fig. HJd). The misfit between the theoretical 
histograms and the observed ones becomes even clearer 
when we consider that all curves in Fig. |5] are independent 
of the position on the solar disk (as they should, because 
they correspond to an isotropic distribution; Eq. [T5)) . while 
the observed ones (Fig. |4]) do change. Choosing a different 
f{B) function in Eq. [T^] does not alter this because the 
histograms of the Stokes profiles would still be independent 
of the position on the solar disk. However, it is possible to 
obtain histograms that vary with the position on the solar 
disk by postulating that the mean magnetic field depends 
upon Bo{Q) or Bo{A). In this case, the different curves 
in Fig. [5] will certainly change with the position on the 
solar disk, yet, in doing so we would implicitly introduce a 
dependence with Q and/or A in the original distribution, 
given by Eq. 1131 and thus it would no longer be isotropic. 
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4.2. Triple-Gaussian distribution functions 

We now consider a distribution function in the local refer- 
ence frame on the solar surface, {gx, ey, Gz}, in which each 
component of the magnetic field is statistically independent 
of the rest and shows a Gaussian-like dependence: 



Pi(B)dB 



dS^dSydS, 



exp 



• exp ■ 



■ exp 



B' 



(15) 



Again, with this definition, the e^-axis is perpendicular 
to the solar surface. The distribution is of course normalized 
to one: 



Vi{B^, By, B,)dB^dBydB, = 1 . (16) 



In the previous equation, B^q, Byo, and B^o refer to 
the mean of the absolute value for each component of the 
magnetic field: 



With this theoretical distribution function of the mag- 
netic field vector we can again (see Section l4?lT) synthesize 
the Stokes profiles, add noise, and finally construct the 
theoretical histograms of the resulting Stokes profiles. 
Equations [TH] and [TO] now entail an explicit dependence on 
Q, unlike the case of an isotropic distribution (Eq. [T5|) . 
Interestingly, if -Bxo — Byo — B^a, the distribution becomes 
isotropic, because Eq. [K] can be re-written in terms of 
the module of the magnetic field vector B (see Eq. (5]) . In 



the following we employ B- 



xO 



B. 



yo 



22G: B^ 



5 G, 



and therefore the resulting distribution is non-isotropic. 
However, as indicated above, we impose B^q = ByQ so 
that the distribution is symmetric in the XY-plane. The 
distributions for B and 7 corresponding to the aforemen- 
tioned values are presented in Figures [T]: and [TJi for two 
different heliocentric angles: 9 = 0° (red curves), 9 = 45° 
(blue). The reason for employing B^o < B^q, ByQ is because 
in the isotropic case (see Fig. |B^) the value of Bq that 
produced a reasonable fit to the linear polarization at disk 
center yielded too much circular polarization and thus, to 
simultaneously fit both, we need to decrease the vertical 
component of the magnetic field vector while keeping the 
horizontal component at the appropriate level. 



Bxo — 



\Bx\Vi{Bx, By, B.yiB^ABydB, 



(17) 



Similar equations can be written for By^ and -B^o- As 
in the previous section, we perform a rotation of angle 9 
around the Gy-axis to transform Eq. [15] into the observer's 
reference frame: {ey.,e.y,e^}. To guarantee that the result 
of this transformation does not depend on the choice of 
rotation axis in the XK-plane we must ensure that i?xo — 
Byo such that Eq. [12] can be written in terms of iJ^ + By. 
With this, we can now rewrite the theoretical distribution 
function as 



dB,dB„dB, 
7'i(B)dB- ^ ' 



B 



'2 



exp ■ 



• exp ■ 



exp ■ 



(g^ sin9-F.BxC0s9)^ | 
(B^ COS 9 - sin 9)2 1 



(18) 



And finally, to solve the radiative transfer equation 
we need to express Equation [TH] in spherical coordinates 
(Eq. [TT|) . This transformation yields 



Pi(B,7,0)dBd7d0 = 



B"^ sin 7di3d7d0 



• exp ■ 



• exp 



• exp ■ 



TT'^iJxO^yoSzO 

[B sin 7 sin 0)^ 1 
J 

(_B cos 7 sin 9 + i? sin 7 cost/) cos 9) ^ 
(_B cos 7 cos 9 — S sin 7 cos sin 9)2 

-^bIq 



(19) 



The first feature to notice is that the distribution of B 
(Fig. [7]:) is independent of the heliocentric angle 9. This 
was to be expected because the module of the magnetic 
field vector B is invariant with respect to rotations. The 
second feature to notice is that the distribution of 7 shows 
a peak at 7 = 90° at disk center (red curve). This is a 
consequence of having imposed < ^^xOi^yO- This peak 
smoothes out as we move toward the poles (blue curves). 



The theoretical histograms for the Stokes profiles 
resulting from the previous distribution (Eq. llOp are 
displayed in Figures [7^ and [71d. As mentioned above, the 
values of -BxOi ^yOi and i?zo were selected such that, for 
signals above 2 x 10~^, the theoretical histograms at disk 
center (9 = 0°; red curves) are comparable to the observed 
ones (Fig. [3]). However, at larger heliocentric angles 
(9 = 45°; blue curves), the mismatch between theoretical 
and observed histograms is very clear, in particular for the 
circular polarization (dashed lines). Our theoretical distri- 
bution function (Eq. 1191) clearly predicts that the circular 
polarization should increase as 9 increases (blue dashed 
lines in Fig. [7]). This happens as a consequence of having 
imposed BxO > ^zO; which means that away from disk 
center, the component of the magnetic field vector that is 
aligned with the observer's line-of-sight, B^ — _B|| , increases 
due to the contribution from i?x (see Eq. [9]). Because the 
slope of the V — B\\ curve is so steep (see Fig. [S]), a small 

increase in i?^ = _B|| translates into a large increase in 
the Stokes V signal. Therefore, the generated amount of 
circular polarization is much larger at larger heliocentric 
angles (dashed-blue curves; 9 = 45°) than at disk center 
(dashed- red curves; 9 = 0°). Likewise, since i?zo < ^xo, 
the component of the magnetic field that is perpendicular 
to the observer's line-of-sight _B_l decreases as 9 increases. 
In this case, however, the linear polarization does not 
decrease as much as the circular polarization increases. 
This is due to the gentler slope of the Q — B± curve (Fig. [5]). 
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Fig. 7. Same as Figure IH] but for the triple- Gaussian theoretical distribution given by Eq. [THl In this case the mean 
value of the three components of the magnetic field vector are different: B^o = Syo = 22 G, and B^o = 5 G. Red lines 
correspond to 8 = 0° (disk center) while green and blue lines correspond to Q = 45°. Solid and dashed lines indicate 
the circular (Stokes V) and linear (Stokes Q and U) polarization. 



Interestingly, Figs|3^-|3lD show that the observed amount 
of circular polarization does not increase, as required by 
the change in the viewing angle, for larger heliocentric 
angles (i.e. as Q increases) in accordance to the observed 
drop in linear polarization. We conclude therefore that the 
differences in the histograms of the observed Stokes profiles 
at different positions on the solar disk cannot be produced 
by the change in the viewing angle Q. This problem is 
not unique to the theoretical distribution given by Eq. [121 
but it will indeed affect any theoretical distribution that 
is prescribed in the local reference frame and features 
Bzo < BxQjByQ. Prescribing a theoretical distribution 
in the local reference frame means that the underlying 
distribution is always the same regardless of the position 
on the solar disk. The only reason that this distribution 
changes is because the angle Q, between the observer's line 
of sight and the vector perpendicular to the solar surface, 
varies with the position on the solar disk. 



4.3. Other distribution functions 



To avoid the problem described in Sect. 14.21 we now 
prescribe a theoretical distribution of the magnetic field 
vector that, already in the local reference frame, depends 
on the latitude A. This will imply, unlike the previous 
distributions in Sections 14.11 and 14.21 that the underlying 
distribution (i.e. in the local reference frame) is different 
at different positions on the solar disk. In particular, the 
A dependence means that the distribution function varies 
toward the poles, but not toward the limbs. The expression 
chosen in this section is 
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• [{Bx cos A + B^ sin A)^ + B^] i [B^ cos B^ sin A]^ 



• exp ■ 



• exp ■ 



\Bz cos A — Bx sin A| 
I [{Ba, cos A + sin A)2 + 3 1 (1 + sin A) 



(20) 



This theoretical distribution is normalized to one as 
it verifies Equation [ini As in the two previous sections, 
we now perform a rotation of angle Q around the ey-axis 
to express Equation [201 hi the observer's reference frame, 
and then we transform into spherical coordinates in the 
observer's reference frame. This is done by applying 
Equations [3] and [11] A close inspection (e.g make for 
instance A = 0) of Eg . [20l shows that again the result does 
not depend on the choice of rotation axis in the A"y-plane. 
Although trivial to obtain, the resulting expression for 
the distribution function in spherical coordinates in the 
observer's reference frame is too long to be written here 
(we provide it for a simplified case below). It is noteworthy 
to mention, however, that even though Eq. [20] depends 
on A, once this equation is expressed in the observer's 
reference frame, it depends both on A and 9: 7^i(B; A, 8). 

The theoretical distributions functions of B and 7 
resulting from Equation [2D] are shown in Figures [S]: and 
[8ji, respectively. Here we considered the following values: 
( = 0.3, = 3, and = 7. Unlike Sections [H] and [121 

and do not correspond to the mean values of B\\ and 
B±. For this reason we refer to them as fi\\ and /3_l instead 
of i3||o and B±^q. In this new distribution, the module 
of the magnetic field vector B is different at different 
latitudes, with a slightly higher mean value at disk center 
(red line; A = 0°) than closer to the poles (green and 
blue lines; A = 30 — 40°). This is a direct consequence of 
prescribing a distribution of the magnetic field vector in 
the local reference frame that depends on A. Otherwise, 
as it occurred in Sects. 14.11 and 14.21 the distribution of B 
would be the same, since the module of the magnetic field 
vector is invariant with respect to rotations. 

We now use the distribution given by Equation [20] 
to solve the radiative transfer equation and produce 
theoretical histograms of the Stokes profiles. These are 
displayed in Figures. [8^ and [SJd. The shape of the curves 
for the circular (dashed) and linear (solid) polarization in 
Figure[Hti are very similar to the observed ones (Fig. [3^). Of 
particular interest is the fact that, in agreement with the 
observations, the histograms of the circular polarization 
are very similar in the three considered positions on the 
solar disk, and feature a peak at signals « 3 — 4 x 10"'^. It 
is also important to mention that the linear polarization 
decreases from A = 0° towards A = 30 — 40°, which also 
agrees with the observations. 

To understand how Equation [20] can fit the observed 
histograms of the Stokes profiles at different latitudes, it is 
convenient to re- write this theoretical distribution function 
assuming that the observer's line-of-sight is along the zero 



meridian. In this case the rotation angle 8 is equal to A, 
thereby simplifying the probability distribution function 
from Eq. [50] into (in spherical coordinates in the observer's 
reference frame) 



i?^ sin'' 7 cos'^ 7 



Vi{B,-i,(t))dBA-id(t> = 



■ exp 



1927r/3|/35 
Bsin7|(l + CsinA) 



(1 + CsinA)^ 



(21) 



,expi-^^|di3d7d0. 

In this Equation [5T] the exponential term that refers to 
B cos 7 (i.e. component of the magnetic field vector that 
is parallel to the observer's line-of-sight) does not depend 
on A and therefore the amount of circular polarization 
generated by this distribution will not change with latitude 
(see dashed lines in Figs. [8^ and [Hd). Moreover, the 
exponential term that contains B sin 7 (i.e. component 
of the magnetic field vector that is perpendicular to the 
observer's line-of-sight) decreases as A increases and thus, 
this distribution function produces less linear polarization 
at higher latitudes (see solid lines in Figs. [5^ and[SjD). 

Despite all these similarities between the theoretical 
and observed distribution of Stokes profiles there are 
still some significant differences, namely: a) there is a 
clear deficit of pixels with sufficient amount of linear 
polarization in the range of signals 2 — 4 x 10"^ 
(Fig. W) or S/R > 5 (Fig. [SJ)) at all latitudes; and b) 
the theoretical histograms show a peak in the circular 
polarization at the 3cr-level at high latitudes (dashed-blue 
and dashed-green curves in Fig. [B^), while the observed 
histograms show no peak in the circular polarization at the 
3(7- level at any disk position (dashed curves in Fig. [3^). 
Overall, however, the fit between the observed and the 
theoretical histograms of the Stokes profiles is clearly 
better for this distribution (Eq. [20l) than for the case of an 
isotropic distribution (Eq.[T21) or a triple Gaussian (Eq.[Tr 



5. Assumptions and limitations 

In the previous subsections we have focused on the effect 
that the probability distribution function of the magnetic 
field vector Vi(B,a) has on the observed histograms 
of the Stokes profiles at different positions on the solar 
disk. In our theoretical analysis a number of simplifying 
assumptions were made to keep the problem tractable. 
Although they have already been pointed out in Section [3] 
we summarize them here to briefly discuss their implica- 
tions. 



— We have assumed that the thermodynamic and mag- 
netic parameters are statistically independent of each 
other. This allowed us to write the total probability 
distribution function in Eq. [H as the product of two 
distinct probability distribution functions. However, as 
dictated by the Lorentz-force term in the momentum 
equation in magnetohydrodynamics, the magnetic 
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Fig. 8. Same as Figure IH] but for the theoretical distribution given by Eq. [501 The parameters employed are /3|| = 3 G 
and /3_L = 7 G. Red lines eorrespond to A = 0° (disk center), while green and blue lines correspond to A = 30° and 
A = 40°, respectively. Solid and dashed lines indicate the circular (Stokes V) and linear (Stokes Q and U) polarization. 



field affects the thermodynamic structure of the solar 
atmosphere, ft is therefore clear that this assumption 
does not fully hold in the solar atmosphere. For 
instance, if we take the commonly accepted picture 
of intergranular lanes harboring more vertical and 
stronger magnetic fields than the granular cells, and 
we consider that intergranular cells have a smoother 
variation of the temperature with optical depth (see i.e. 
Fig. 3 in Borrero & Bellot Rubio 2002), we could then 
have postulated a correlation between the magnetic 
field vector B and the gradient of the source function 
with optical depth 5*1, which is contained in T (Eq. [3]). 
Indeed, the higher Si is, the stronger will be the 
polarization profiles Q, U, and V (see Eq. 9.45 in del 
Toro Iniesta 2003). These correlations could potentially 
cause the observed histograms of the Stokes profiles 
(Figure [3]) to vary with the heliocentric angle, even if 
the underlying distribution of the magnetic field vector 
does not depend on Q. Therefore it is important to 
investigate what an effect they have before conclu- 
sively proving that the distribution of the magnetic 
field vector is not isotropic (Sect. I4.ip . or that the 
differences in the observed histograms of the Stokes 



profiles are not due to the viewing angle (Sect. 14. 2|) . 
Unfortunately, the aforementioned correlations are not 
known for the solar internetwork simply because it is 
not clear how magnetic fields are distributed here. In 
the future we will explore this question by employing 
3D numerical simulations of the solar atmosphere, 
because they provide correlations between B and T 
that are compatible with the MHD equations. 

— We have also assumed that the probability distribu- 
tion function of the thermodynamic and kinematic 
parameters, 7^2(T,Vios), does not depend on the 
position on the solar disk. Kinematic parameters (i.e. 
line-of-sight velocity Vios) do not influence our study, 
since they have no effect on the amplitude of the 
Stokes profiles in Fig. S) The same can be argued 
about other thermodynamic parameters in T, such 
as the source function at the observer's 5*0 (affects 
only Stokes /), and damping parameter a (affects 
mostly the line width but not its amplitude). By 
far, the most important thermodynamic parameters 
affecting the amplitude of the Stokes profiles under 
the Milne-Eddington approximation are the gradient 



11 



J.M. Borrero and P. Kobel: Inferring the magnetic field vector in the quiet Sun 



of the source function with optical depth 5*1 and the 
continuum-to-line-center absorption coefficient 770. In 
a ID atmosphere both these parameters are known 
to decrease as 8 increases, because the line-of-sight 
samples a thinner vertical-portion of the atmosphere. 
However, since the dependence of the polarization 
profiles with Si and 770 are identical (see Eqs. 8.14, 
8.15 and 9.44 in del Toro Iniesta 2003), one would 
expect the same drop with increasing heliocentric angle 
angle in the amplitude of the circular polarization 
profiles (Stokes V) and linear polarization profiles 
(Stokes Q and U). However, Figure U shows that the 
linear and circular polarization profiles (solid-color and 
dashed-color lines) behave differently, and therefore 
we can rule out the variations of 5*1 and/or 770 with 
B as being responsible for the observed histograms in 
the Stokes profiles. Of course, this would change in a 
3D atmosphere, where the line-of-sight pierces through 
different inhomogeneous atmospheric layers, thereby 
opening the door for the possibility of S'l and or 770 
to affect the linear and circular polarization profiles 
differentljQ. 

Adopting a Milne-Eddington atmosphere also implies 
that we are assuming that the magnetic field vector B 
does not vary with optical depth Tc in the photosphere. 
This can have important consequences, since at larger 
heliocentric angles the spectral line samples higher 
atmospheric layers than at disk center, where the 
probability distribution function of the magnetic field 
vector can be different. Employing the widely used 
ID HOLMU model (Holweger & MiiUer 1974) we 
calculated that the continuum level Tc = 1 rises by 
approximately 20 Km from disk center = 0° (map 
A; Sect.mil) to 9 = 45° (maps B and C; Sects. and 
12. 3p . Since this vertical shift of the continuum level is 
rather small, we could argue that the histograms of the 
Stokes profiles in Fig. |3]are not affected by this effect. 
However, the value of 20 km should be considered 
only as a lower limit since a ID model does not take 
into account the horizontal inhomogeneities present in 
the solar atmosphere. To properly account for this ef- 
fect, more sophisticated 3D models should be employed. 

Finally, we have considered a = 1 in our analysis (see 
Sect. EJ. This is equivalent to considering that, at 
the resolution of the Hinode/SP instrument (0.32"; 
Sect, the magnetic structures are spatially resolved. 
This is, of course, highly unlikely, and therefore it 
would be important to drop this assumption in the 
future. Its importance can only be quantified with 
additional assumptions about the scale-distribution of 
the magnetic structures in the solar photosphere. This 
topic is, in itself, as controversial as the distribution of 
the magnetic field strength and inclination, which is 
the reason why we have refrained from addressing it 
here. Although employing 3D MHD simulations would 
certainly help to drop the a = 1 assumption, we are 
cautious about it since it is not clear whether these 



simulations are reliable at the smallest physical scales 
(Sanchez Almeida 2006). 



6. Discussion and conclusions 

The histograms of the observed Stokes profiles at different 
positions on the solar disk (Fig. 2]) are clearly different 
from each other. One possible interpretation for this 
is that the distribution of the magnetic field vector in 
the solar internetwork is not isotropic. We explored this 
possibility in Section 14.11 where we employed an isotropic 
probability distribution of the magnetic field vector. This 
distribution yielded, as expected, the same distribution of 
Stokes profiles at all positions on the solar disk (Fig. 
Martinez Gonzalez et al. (2008) have also presented similar 
histograms but employing the Stokes profiles from the Fe I 
line pair at 1.56 /im (observed with the TIP2 instrument; 
Martinez Fillet et al. 1999). Their histograms (see their 
Figure 2) showed no clear variation with the heliocentric 
angle, which lead them to conclude that the distribution 
of the magnetic field vector in the quiet Sun was isotropic. 
Interestingly, these authors also mentioned after a more 
detailed analysis that there could indeed be a dependence 
of the histograms with the heliocentric angle (as indeed we 
find here). 

In addition to Martmez Gonzalez et al. (2008), a 
number of works have also argued in favor of an isotropic 
distribution of magnetic fields in the internetwork. In 
particular, Asensio Ramos (2009) and Stenflo (2010), 
employing two different approaches, both concluded that 
for very weak magnetic fields {B — >• 0) the distribution 
becomes isotropic. With our present data we cannot argue 
against or in favor of this interpretation. The main reason 
for this is that, as discussed in Section [31 any distribution 
for the magnetic field vector where B± has a peak below 
40-70 will produce linear polarization profiles that are 
dominated by noise (3(T-level or S/R = 3). Therefore 
our current approach (described in Section 3]) cannot be 
employed to discern the underlying distribution of the 
magnetic field vector from these profiles dominated by 
noise. However, it can be employed to establish that the 
number of pixels that would follow this hypothetically 
isotropic distribution cannot be much larger than 30 % of 
the pixels in the internetwork, since this is the amount of 
pixels that show a peak at the 3<7-level in the polarization 
profiles (see Fig. 2^). For signals above > 2 x 10~^, the 
histograms of the Stokes profiles deviate significantly from 
the ones predicted by an isotropic distribution, and thus 
we can establish that here the distribution of the magnetic 
field vector cannot be isotropic. 

We can use a different argument to further clarify the 
previous point. Our theoretical distributions in Section 2] 
apply to all possible values of the module of the magnetic 
field vector. However, we could have employed distributions 
pieced together in the following form: 



^ Outside the Milne-Eddington approximation, a similar argu- 
ment could be made in terms of the stratification with optical 
depth of the temperature in the solar photosphere Tijc), instead 
of Si and 770. 



7'i(B)dB 



7'a(B)dB, \iB<B* 
n(B)dB, \iB>B*, 



(22) 



where Va could hypothetically correspond to an isotropic 
distribution for weak fields: B < B* . This would explain 
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the 3CT-peak in the linear polarization in Figure 0^ (dashed 
lines). In addition, Vb could be a distribution, valid for 
larger fields B > B*, that would fit the tails of the 
histogram. The distribution given by Equation [52] does not 
need to be discontinuous because it could be prescribed 
such that Va{B*) = VbiB*). 

In Section 14.21 we employed a triple Gaussian (one for 
each component of the magnetic field vector) distribution 
function and found that, under this assumption and at 
disk center, the best fit to the observed histograms of the 
Stokes profiles is produced by a distribution in which the 
mean value of the magnetic field vector component that is 
parallel to the solar surface is lower than the mean value of 
the magnetic field vector component that is perpendicular 

to the solar surface: Bzq < sj S^o + -^yo • This yields a 
distribution function where the magnetic field vector is 
highly inclined, in agreement with previous findings from 
Orozco et al. (2007a, 2007b) and Lites et al. (2007, 2008). 
However, this distribution does not fit well the histograms 
of the Stokes profiles at other positions on the solar disk. 
In fact, in that section we found that it is not possible to fit 
the observed histograms for the Stokes profiles at different 
heliocentric angles employing a theoretical distribution 
function for the magnetic field vector prescribed in the 
local reference frame that only changes due to the viewing 
angle Q. The reason for this is that, for an underlying 
distribution where the magnetic field vector is mostly 
horizontal at = 0° (disk center), the amount of linear 
polarization slightly decreases when O increases, while 
the amount of circular polarization would significantly 
increase as <d increases. However, the observed histograms 
of the Stokes profiles (Fig. 0]) show that, although the 
amount of linear polarization decreases when Q increases, 
the circular polarization does not particularly increase 
(see also discussion in Lites et al. 2008). This cannot be 
explained in terms of a simple rotation of the viewing angle 
6, and therefore we interpreted this fact, in Section as 
proof that the underlying (i.e. in the local reference frame) 
distribution of the magnetic field vector must depend on 
the position on the solar disk. 

Under the assumption that the distribution of the 
underlying magnetic field vector depends on the latitude 
A (see Sect. 14. 3L we were able to find a theoretical 
distribution of the magnetic field vector (Eq. [20)1 that fits 
quite well the observed histograms of the Stokes profiles 
at different positions on the solar disk (Figure [8]). Among 
other properties, this distribution features a magnetic 
field whose mean value decreases toward the poles. We 
note here that this does not mean that this is the real 
distribution for the magnetic field vector present in the 
quiet Sun. One reason for this is that the fit is far from 
perfect (see discrepancies mentioned in Sect. l4.3"| . but most 
importantly, that we do not know whether this solution is 
unique because there can be other theoretical distributions 
that fit the observed Stokes profiles equally well, or even 
better. More work is indeed needed to confirm or rule 
out Eq. [50] as the real distribution of the magnetic field 
vector present in the Sun. In particular, a better fit to the 
observed histograms of the Stokes profiles is desirable. In 
addition, it is important to have maps at more latitudes to 



further constrain the possible distribution functions. 

Moreover, it is important to bear in mind that the 
conclusions above are not necessarily the only possible 
interpretations, because postulating a probability dis- 
tribution function of the thermodynamic and kinematic 
parameters, ■P2(T, Mos), that varies with the heliocentric 
angle 0, or postulating a correlation between the thermo- 
dynamic (T) and magnetic (B) parameters might also help 
explain the observed differences between the histograms 
of the Stokes profiles at different positions on the solar 
disk. Another effect that has not been accounted for is 
that magnetic field vector can vary with optical depth in 
the solar photosphere. Since the Stokes profiles sample 
increasingly higher atmospheric layers as the heliocentric 
angle increases, the distribution of the magnetic field 
vector can be different for different values of 0, even if 
the probability distribution of the magnetic field vector 
is the same at all positions on the solar disk at a fixed 
geometrical depth. All these effects could be properly 
accounted for by means of 3D MHD simulations of the 
solar photosphere (Schiissler & Vogler 2008; Steiner et al. 
2008, 2009; Danilovic et al. 2010). 

In the future we expect to employ such simulations 
to either rule out or confirm our results in this paper. 
Consequently, our conclusions at this point should be 
regarded as preliminary only. Instead, the main purpose 
in this paper is to illustrate the methodology detailed in 
Sections [3] and [3] to study the distribution of the magnetic 
field vector in the quiet Sun, by directly inverting the 
histograms of the Stokes profiles in entire maps instead of 
inverting the Stokes profiles at each spatial position in a 
given map. Our method has great potential to investigate 
several aspects of the photospheric magnetism in the solar 
internetwork. For instance, it can be used, as in Sect. 14.31 
to confirm whether the mean value of the distribution of 
the magnetic field vector changes from disk center toward 
the poles (cf. Zwang 1987; Ito et al. 2010). This will have 
important consequences for theoretical models that explain 
the torsional oscillations in the butterfiy diagram in terms 
of a geostrophic flow model (Spruit 2003), which requires 
a significant amount of magnetic flux at high latitudes 
at the beginning of the sunspot cycle. In addition, and 
although in this work we have restricted ourselves to 
variations in latitude (A), additional observations from 
disk center toward the solar limbs could be employed 
to investigate whether the properties of the magnetic 
field in the internetwork change also in longitude. This is 
already predicted by non-axisymmetric dynamo models 
(Moss 1991, Moss et al. 1999, Bigazzi & Ruzmaikin 2004, 
Charbonneau 2005) and can provide important clues about 
the strength of the differential rotation (Riidiger & Elstner 
1994; Zhang et al. 2003). 
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